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Abstract 

Quantum parameter estimation, the ability to precisely obtain a classical value in a quantum system, is very 
important to many key quantum technologies. Many of these technologies rely on an optical probe, either coherent 
or squeezed states to make a precise measurement of a parameter ultimately limited by quantum mechanics. We 
use this technique to theoretically model, simulate and validate by experiment the measurement and precise 
estimation of the position of a cavity mirror. In non-resonant systems, the achieved estimation enhancement from 
quantum smoothing over optimal filtering has not exceeded a factor two, even when squeezed state probes were 
used. Using a coherent state probe, we show that using quantum smoothing on a mechanically resonant structure 
driven by a resonant forcing function can result significantly greater improvement in parameter estimation than 
with non-resonant systems. In this work, we show that it is possible to achieve a smoothing improvement by 
a factor in excess of three times over optimal filtering. By using intra-cavity light as the probe we obtain finer 
precision than has been achieved with the equivalent quantum resources in free-space. 

PACS numbers: 42.50.Dv; 03.65.Ta: 03.67.-a 
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Background 

Introduction 

The field of quantum metrology can be described as using quantum resources to enhance measurement 
precision beyond that achievable with purely classical resources. There are a number of resources that 
are available such as entanglement [T], superposition [2] and squeezing [3]. There are also tools such as 
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adaptive feedback [4] and quantum smoothing [5] to further exploit the quantum enhancement. Quantum 
parameter estimation (QPE) is a related discipline that is focussed more specifically on precisely estimating 
the classical parameters of a quantum system. The importance of QPE to fields such as gravitational wave 
detection [6], quantum metrology laiH], quantum control [9] and opto-mechanical force unmi] sensing has 
been well established. Technological evolution in recent times has seen an increase in the range of pertinent 
architectures where quantum mechanical effects have become relevant [niiiain]. There have also been 
experimental demonstrations of key advances in QPE. Eor example, in optical phase estimation we saw 
successive lowering of achievable mean square estimation error by the use of adaptive feedback m and 
adaptive feedback was combined with smoothing m to achieve a further reduction. With the addition of 
phase quadrature squeezing the limit was once more lowered m- A recent extension of these QPE techniques 
to a more macroscopic domain uses an optical probe beam to obtain an estimate of the position, momentum 
and force acting on a free-space mirror m- 

As was shown in [5], an increase in estimation precision relative to filtering is expected when quantum 
smoothing is used. Previous work here has only considered first-order forcing noise processes with non¬ 
resonant interactions between the forcing functions and the system. Eor such set-ups, results to date have 
yet to show a greater than two improvement of the smoothed estimate over the filtered equivalent. An 
interesting open question therefore remains as to whether this factor of two improvement is an upper limit 
for more complicated systems. So in this work we consider a higher order forcing function that is Lorentzian 
in frequency. Additionally, we consider resonant interactions between the forcing function and the system 
(mirror) with the centre frequency of the Lorentzian aligned with the peak of a mechanical resonance. Here 
our theory suggests that for more realistic resonant systems driven by less benign processes, the factor of 
two improvement with smoothing can be improved on significantly. We present theory and simulations 
results showing a greater than two smoothing improvement over the equivalent optimal filtered estimate 
obtained. The results of the simulations both verify and extend beyond the theoretical analysis and we 
present experimental results to verify the simulations. 


Theory - optics 

To date the experimental demonstrations of smoothing have focussed on systems where the probe beam 
(even when quantum enhanced) interactions are in free-space. It is relatively well known that optical cav¬ 
ities can be used to enhance measurement precision. In the context of this work the strong intra-cavity 
field in an optical cavity provides more photon interaction with the parameter to be estimated. As each 
photon potentially probes the parameter many times the cumulative effect gives higher sensitivity with¬ 
out need for extra photon resources. Because the experimental validation makes use of the enhancements 
in sensitivity achievable by the use of optical cavities, the theory and simulation assume an intra-cavity probe. 


Using the formalisms in [18], a single ended cavity, see figure is described 


in terms of optical fields as 


d — —K,a H- iAa + \/2^Ain -h 


( 1 ) 


where Ain is the input field and a is the cavity field. Here k, = K,a i^a is the half width half maximum 

(HWHM) cavity decay rate and Kia is the intra-cavity loss term. We use the standard approach of separating 
AC and DC terms, i.e. A = a SA and note that the loss term (Ai) has no DC component. After 
solving Q and its conjugate for steady state (d = 0 & d* = 0) and applying the boundary condition 
(<^out = — ain) for the output coupler we obtain 
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We use the standard quadrature definitions 


i^lut + ^out) and - Aout) and assume 
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Figure 1: Single ended cavity model used in derivations. 


that Ain is real. We apply an AC forcing function to a cavity mirror via a piezoelectric transducer (PZT) 
(bottom right figure that varies the mirror position. The goal is to estimate that forcing function with the 
smallest possible mean square error (MSE). The applied signal varies (or detunes) the resonant frequency of 
the cavity about an average value set by cavity length. This is manifested in equation 0 via the non-linear 
cavity detuning term A, which is zero when the optical frequency is equal to the cavity’s resonant frequency. 
The result is a non zero detuning term and hence a non zero signal on the phase quadrature at the applied 
AC frequency. To account for this, we also separate the non-linear cavity detuning term into average and 
fluctuating terms giving A = A + C(t), where C(t) accounts for our applied AC forcing function (see equations 
0 and §)• The DC component of the detuning term A is assumed to be zero, meaning the DC component 
of is neglected and the DC field of interest is 


x+ ^ 



( 4 ) 


Now we address the fluctuating terms in Q. We cannot assume steady state so we move to the Fourier 
domain to solve the differential equation. Using the relation that = ioox and substituting the DC 

solutions for a and a* as necessary we obtain 


iujbA — — Kj5A T id^cx -h Ain T (5) 

-h y/2i^5A\^ + ^/2ki5A\. (6) 


As we are interested in low frequencies, we assume that uj k. After applying the boundary conditions and 
some relatively straight forward algebra, the output quadratures are 
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The fluctuating component of the detuning term (C) includes the PZT response to higher frequency pertur¬ 
bations, i.e. the applied forcing function. It therefore makes sense to think about the output field, which is 
our probe of the mirror position, from a signal and noise sense. We assume that only the phase quadrature 
of the probe is measured. The phase quadrature (equation 0 ) can been seen to consist of components that 
are quantum fluctuations (first two terms) and a component that is a function of the detuning parameter 
(last term). The detuning (C(<^)) of fho optical resonant frequency of the cavity is a function of the cavity 
length changing with mirror displacement. The mirror is coupled to the PZT and so the magnitude of the 
displacement depends on the PZT’s frequency response. At a mechanical resonance of the PZT, a greater 
displacement will be imparted on the mirror for a given forcing function. Therefore ("(u;) is not constant 
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with frequency but varies as a function of the PZT transfer function. By contrast {uj) and {uj) 

represent quantum vacuum fluctuations which are constant with frequency. 


Theory - smoother 

In this subsection, we develop the theory to predict how much improvement can be expected by using 
smoothing to obtain our estimate compared to filtering. In the process, we will derive expressions for the 
optimal smoothed MSE and the transfer function of the optimal smoother. We take a block diagram approach 
and consider the optical part of the system as a generic plant with input and output signals (I/O), as shown 
in figure At this stage it is more intuitive to consider signal voltages rather than perturbations in metres. 
The system definition will be recast later to derive the optimal smoother for estimating the mirror position 
in metres. From figure we define our system as 



Figure 2: The block diagram used for the derivation of the controller and the voltage estimator. 


Vy{t) 




( 9 ) 

( 10 ) 


where v^{t) is the noiseless (unmeasurable) signal due to the plant disturbance, u^(t) is measurement noise 
(in our case dominated by quantum noise) and Vy{t) is the measured output of the homodyne detector. Also 
h{t) is the transfer function of the entire optical system, Vf{t) is the applied forcing function, Vc{t) is the 
control signal and * represents a convolution operation. We now develop the optimal smoother and establish 
the MSE of the smoothed estimate for the system defined by and ( prl ). Assuming that we can obtain a 
stable controller hc{t): 


Vc{t) = hc{t)^Vy{t) 

v^{t) = h{t) * (u/(t) - hc{t) * [Vy,{t) + 


We continue in the Fourier frequency domain where equation ([T^ becomes 

h{uj)vf{uj) — h{uj)hc{oo)vy{uj) 


,{uj) = 


1 T h(cij)hc(cij) 


( 11 ) 

( 12 ) 


(13) 


As there may be uncertainty in the system parameters, we prefer the smoothed estimate {vs{t)) to be 
minimally impacted by a sub-optimal control signal (uc(t)). This controller independence is achieved by the 
addition of the h{t) block in the estimator box of figure so that 

Vz{t) = h{t) ^ Vf{t) — h{t) ^ Vc{t)Vy{t) ^ h{t) ^ Vc{t). (14) 
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The Vc{t) terms in (14) cancel giving 


Vz{t) = h(t)^Vf{t) +Vr^{t). (15) 

Here we note that the optimal estimation theory assumes a known transfer function h{uj)^ which we measured 
experimentally. Although this assumption suffices for our purpose, any mismatch between our assumed model 
and reality will lead to an increase in the actual error. Many techniques are available to address this potential 
problem m- We now derive the transfer function for the optimal smoother hs{oj) to estimate Vf{uj), 

Vs{uj) = hs{uj)vz{uj) (16) 


and define the estimation error as 


Av(cj) = Vs(cj) — Vf(cj) = Vf(cj)[hs(cj)h(cj) — 1] + hs(cj)vy(cj). (17) 

The power spectral density of the error signal is 

S'a«(w) = lhs(oj)h(cj) - H^Svj:(oj) + hs(oj)I^Sv^(oj). (18) 

The mean square error (MSE) is defined in the normal way and can be expressed in the frequency domain 
using Parseval’s theorem [20] 


/ OO 7 

4^S'a«(w). 
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From equation (19), the mean square error is 


(19) 


( 20 ) 


We minimise the MSE by finding the functional derivative of equation (20) with respect to hs{u)) to obtain 
the optimal smoother 


ha{uj) = 


h*{u})Sv„ (w) 


\h{u})\^Svj{uj) + Sv^{uj)' 


By substituting (21) into (20), the mean square error for the optimal estimate of Vf{t) is found to be 

du 5„,(w)5„,(w) 




-OO 27r |/i(w)|2S'„.(w) + Sv (u) 


( 21 ) 


( 22 ) 


To recast the smoother to estimate mirror position we insert a conversion block and redefine the transfer 
function, see figure It is then just a matter of reworking the above derivation with the new system. After 
reworking the algebra, we find that the optimal smoother for position estimation is 


where Sxj = Svf ■ and h'{Lo) 

found to be 


\h'{L0)\^S,^+S,^ 


(23) 


= h{uj)/Ap^j.. Prom which the optimal mean square position error is 


du) SxfSi)^ 

J-oo 27r \h'{u})\‘^Sxf + 


(24) 


This can be shown simply and conveniently to be = APp^^Cy. Here Ap^p is a constant that relates the 
voltage applied to the PZT to the physical mirror displacement and was measured to be Ap^p « 6.3 x 10“^ 
[m/V], 
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Figure 3: The block diagram with additional gain block used for the derivation of the position estimator. 
Here the transfer function h{t) is modified to h'{t) that assumes position, Xe{t) (units of metres), rather 
than voltage, Ue(t), as an input. As such the smoothed output, Xs{t), is similarly a position estimate. 


Plant and forcing function 

So far, no assumptions have been made about the dynamics of the forcing function or the plant. Of interest 
in this work is the MSE when those dynamics are non-trivial. In this subsection we describe both the plant 
and the forcing function used in the theory, simulation and later the experiment. Starting with the plant, 
the model used in the theory and the simulation is based on modelling a true experimental system. The 
transfer function of a true physical system was measured using a dynamic signal analyser. We then modelled 
the dominant resonance, see h{uj) in table The magnitude and phase plots of the measured system and 
the model are shown in figure The plot shows that there is a time delay (identified by the constant phase 
lag super-imposed on the other features in the lower plot of figure |^. This delay is included in the model 
but not compensated for in either the controller or the smoother. It can be seen that only the dominant 
mechanical resonance at cj ^ 27r x 7640 rad/s is accounted for in the model. The other resonances are 
approximately 20 dB down and are taken to be not excited by the forcing function. 

The system is driven via a cavity mirror with a Lorentzian forcing function {vf{t)) as shown at figure 

Again, measurements of a physical system are used as the basis of the model. Figure shows the PSDs 
of both the simulated (blue) and experimental (red) forcing functions. The model used in the controller 
and system design is shown in table {Sf{uj)) and is accurately represented by the blue plot in figure]^ 
The experimental forcing function is also a good match to the theory for the frequency range of interest 
(< IbkHz). It is apparent that the process used in the experiment has higher order (odd) harmonics that 
are not accounted for in the models. These are an artefact of the experimental generation of the forcing 
function. The higher order harmonics do not excite the system because both the plant and the controller 
are heavily attenuated at these frequencies. 


Results and Discussion 

Simulation 

We developed a numerical simulator to test the theory and provide a baseline against which an experimental 
testbed can be compared. The numerical results can also be used to inform subsequent experiments. The 
simulation was done using Simulink and the parameter values are shown in Table The input and mea¬ 
surement noise processes were entered as floating point arrays from the workspace. The plant {h{t)) and 


6 














































Table 1: System parameters for simulation and experimental validation of the simulator. 


Parameter 

Simulation 

Experiment 

Description 

h{uj) 

CiS-\-C2U^rn , —.st 
s2+/3s+a;^ 

CiS-\-C2idm 

S^+pS+UJ^ 

Plant transfer func¬ 
tion 

Sfico) 

Q \ 1 1 1 1 

<3 r 1 1 1 1 

Forcing function PSD 

2 [(a;-c^i)2+72 ' (uj+Wi)^+'r^ \ 


R 

7.7 X 10-11 

7.7 X 10-11 

Measurement noise 
magnitude term 

where R5{t — t') = 
a{r]{t),r]{t)), r]{t) is 
white Gaussian noise 

Q 

7.4 X 10-2 

7.4 X 10-2 

Forcing function 

magnitude term 

where Q5{t — t') = 

£.{t) is 

white Gaussian noise 

7 

1333 

1333 

Forcing function 

damping factor 

LUm 

27r • 7930 

27r • 7930 

Frequency of PZT 

resonance 

iUi 

27r ■ 7930 

27r ■ 7930 

Frequency of forcing 
function resonance 

Cl 

131 

131 

Gonstant 

C2 

196 

196 

Gonstant 

P 

2494 

2494 

PZT resonance 

damping factor 

T 

0 and 18.5 xlO-® 

system 

Time delay 

F 

250 kS/s 

250 kS/s 

Sample rate 

N 

215 

216 

Number of samples 

Averages 

21 

5 

number of data sets 
averaged 


controller {hc{t)) transfer functions were implemented using transfer function blocks with the numerator and 
denominator coefficients extracted from the workspace. The Simulink model provides the parameter Vz{t) 
(see figure for the smoother. In the experimental validation discussed later a cut-down version of the 
Simulink model was used to process the experimental data to obtain Vz{t) from the recorded experimental 
values of Vy{t) and Vc{t). The smoothing (both simulation and experimental) was implemented in the fre¬ 
quency domain using the Fourier transforms of the relevant parameters from the workspace and the Simulink 
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Figure 4: Magnitude and phase plots of the measured system transfer function {h{uj)) and the function fit 
used for the simulation. 


model. The goal of the simulation was to find whether there exists a range of parameters (preferably exper¬ 
imentally feasible) that allow for a greater than two improvement over the optimal filtered estimate. 

The simulation results are shown in figure]^ This figure shows the smoothing improvement factor (S) for 
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Figure 5: Power spectral densities of the experimental and simulated forcing functions plotted against 
frequency. 






























































































































a coherent state as a function of the parameters Q and 7 of the input forcing function {Sf{uj) see Table . 
The smoothing improvement factor is defined as 


S = 


^filt 


(25) 


where is the smoothed MSE (see equation (24)) and Cfnt is the optimal filtered MSE error found numer¬ 
ically using the optimal Kalman-Bucy filter covariance matrix m- The input forcing function parameter 
varied in the upper plot is 7 , which is varied from 500 to 3082 for a fixed Q of 2.35. In the lower plot we vary 
Q from 9.3 x 10“^ to 2.35 for a 7 of 500. The lines are included to guide the eye between the data points. 
The error bars are the standard deviation of 21 separate simulations for each data point. The blue dashed 


line on each plot shows the theory using equation (24) and the equivalent filtered MSE for the respective 


parameters at the data points. The red dash-dot line shows the result for the simulation. The simulated 
data utilised a controller that was designed using the linear quadratic Gaussian (LQG) [ 20 ] methodology 
with fixed central values for the Q (0.57) and 7 (1193). The other parameters used in the LQG controller 
design are shown in Tablewith LQG design parameters /tq = 1 and xq = 2.4. 

It is quite clear in figure that there is good agreement between simulation and theory, with most data 
points agreeing within error bars. It is noted that the model assumes that the optical cavity remains linear 
and as such does not account for large values of the non-linear detuning term (A in equation ^). The 
theory and simulation show that a greater than two smoothing enhancement for a resonant process acting 
on a mechanically resonant structure is achievable for a wide range of parameters. 


Experiment 

The goal of the experimental results presented in this paper are to validate that the theory and simulation 
results appropriately reflect a physically reasonable experimental system. The experimental set up is shown 
in figure We use the 1064 nm output of an Innolight “Diabolo” doubled NdYAG laser as the primary 
optical frequency for the cavity. The 1064 nm beam is spatially filtered using an MCC. After the MGG, 35.6 
mW is split off at a 100:1 ratio with a polarising beam splitter (PBS) for use as the local oscillator (LO) 
for balanced homodyne measurement. The remaining light of approximately 400 /iW is phase modulated at 
199 MHz (REl on figure]^ to create a weak coherent state (n ~ 10^ photons per second). This modulated 
coherent state is used as the input to a single ended bow-tie cavity, with an ESR of 199 MHz [ 21 ]. The 
cavity is locked using dither locking [ 22 I [23] with a frequency of 1.322047 MHz (RE2 on figure applied 
through the EOM and detected at PD3 on figure [^ This dither signal is demodulated and a low frequency 
PI controller is used to maintain the DG frequency locking of the cavity (A in equation 0 ). The cavity 
output (Sig) is sent to a spatial balanced homodyne detector with a fringe visibility of 95.9% and 97.9% for 
each detector respectively (averaged in the quantum efficiency calculation). The low frequency (LE) output 
of the homodyne detector via a PI controller is used to lock the detection to the phase quadrature of the 
signal. The high frequency (HE) homodyne output is demodulated using a 199 MHz (REl) radio frequency 
(RE) LO, an RE mixer (Mixl) and a low pass filter (LPE). This demodulated signal is then feed into the 
feedback filter (EBE). The input signal to the EBE is stored using an Acqiris data acquisition system 

with a sampling rate (E) of 250 kS/s for post processing. The output of the EBE is captured before 

an attenuator (Attn) for better signal to noise ratio and is also stored for post processing. The output of the 
attenuator {vc{t)) is added to the DG lock signal and the applied forcing function (Noise). It is then amplified 
by a high voltage amplifier (HV amp) and applied to the PZT attached to the cavity mirror to be estimated. 
The forcing function signal {Sf{uj) in Tableand also figure]^ is generated by amplitude modulating an 
Ornstein-Uhlenbeck (OU) process generated with an operational amplifier circuit and a white noise generator 
with a carrier of frequency cjc- This signal {vf{t)) is stored for comparison with the estimation and applied 
to the cavity mirror PZT. The signal Vf{t) imposes the variation of mirror position that is to be estimated 
and is the reference for calculation of the MSE. As the data acquisition system has a 50 D input impedance 
all acquired signals are buffered with unity gain operational amplifier circuits so that the acquisition has 
minimal effect on the voltage levels. Additionally anti-aliasing filters were used on all channels of the data 
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Input noise parameter Q 


Figure 6 : Smoothing improvement factor (S) for a coherent state for varying input forcing function {Sf{uj) 
parameters Q and 7 . 


acquisition system and the laser’s resonant relaxation oscillation was suppressed by its noise eater. 

In the previous sub-section, we have presented the results of the theory and simulation. These results 
show the predicted smoothing enhancement consistent with theory and now we use the experimental data as 
validation of the theory. The measured experimental parameters are summarised and compared to those of 
the simulation in Table A controller was constructed using analogue electronics and standard controller 
design techniques [24]. Whilst the design process suggested the controller was stable, the control input was 
also variably attenuated for additional safety. The MSE results presented here (see figure lines included 
to guide the eye) are plotted as a function of this variable attenuation. There are two significant technical 
differences between the experiment and the simulation that need to be considered. The first mentioned 
earlier is the time delay, this is simply accounted for in the Simulink model with a delay block. The second, 
also mentioned earlier (see figure [^, is the higher order harmonics on the captured forcing function. These 
harmonics are below the system noise floor so their effects are negligible in terms of the estimation. As they 
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Figure 7: The full experimental setup for the coherent state cavity mirror position estimation. Note mode 
cleaning cavity (MCC) locking circuitry is omitted. 


are not modelled in either the theory or the simulation the experimental data was corrected as follows and in 
figure]^ for the validation. To obtain the actual MSE from equation (24) the integration limits are infinite. 
In practice this is impossible due to finite sampling rates. So the integration limits were set at ±125 kHz 
limited by our 250 kS/s sampling rate. The contribution of frequencies greater than 125 kHz was found to 
be of the order of 1% for this system. The 125 kHz limit is thus considered reasonable for the theory and 
simulation comparisons, but it includes the harmonics in the experimental case. The experimental MSE is 
artificially inflated relative to the simulation where the forcing function has no harmonics, shown by the 
separation of the experimental (blue triangles) and the simulated (green diamonds) MSEs in figure To 
correct and allow for a fair comparison, we truncate the integration range to ±15 kHz for both. The small 
separation of the MSE (green diamonds) and the reduced frequency MSE (black squares) in figure shows 
that for the simulation, this truncation has only a minor impact. However, due to the removal of the un¬ 
modelled harmonics, the impact of the truncation is much greater in the experimental MSE (pink inverted 
triangles). The reduced frequency MSEs include only modelled data and so are suitable for the experimental 
validation of the simulation. Eigure shows that the corrected MSEs for the simulation (black squares) 
and the experiment (pink inverted triangles) are consistent, thus validating the simulation. The MSEs also 
show a degree of independence from the control input (horizontal with varied control gain) as a result of the 
deliberate cancellation of the control signal dependence in the smoother design. Controller independence is 
likely to be a useful feature in situations where uncertainty in the system model exists and will be further 
investigated in future work. 


Einally we use the experimentally validated simulator to confirm the theory in the current experimental 
parameter regime. Eigure(error bars omitted for ease of viewing) shows the MSE for the simulation (green 
diamonds) as compared to the theoretical smoothed MSE (red solid). The offset can be explained by the 
fact that the smoother does not compensate for the time delay. By removing the time delay block from the 
simulator (cyan circles) this offset is removed and the simulated MSEs become consistent with the smoothed 


11 


























Figure 8: The actual MSEs (blue triangles and green diamonds) and the frequency truncated MSEs (pink 
inverted triangles and black squares) in position for both experiment and simulator as a function of feedback 
gain. 
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Eigure 9: Mean square position error for both experiment and simulator as a function of feedback gain. 


theory. There is a gradual decline in precision for smaller attenuation values for both the experimental 
and the simulation results (see figure]^, this may be due to controller sub-optimally starting to make the 
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system go unstable. The final data in this plot is the MSE optimal filtered estimate (dashed purple) which 
is included for comparison. It is evident that the smoothed estimate is comparable with but not better than 
the system with time delay. However, the no time delay estimate is clearly better than the filtered equivalent 
but has a smoothing improvement factor of less than two but this is expected in this parameter space. One 
final point of interest is in the actual value of the MSE achieved which is approximately 1 x 10“^^ m^ with 
\a\‘^ ^ 10^ s~^. This is comparable with the position MSE achieved in m with squeezing enhancement 
and \a\‘^ ~ 10^ s~^. Whilst this result in itself is not surprising as it is known that optical cavities provide 
additional sensitivity, it is still a good outcome. With future improvements to the system, we expect to see 
further significant lowering of this coherent state MSE. 


Conclusions 

We have developed theory describing resonance enhanced mirror position estimation of a cavity mirror using 
quantum smoothing. This theory has been used to design a numerical simulation model, which we have 
experimentally validated. We have demonstrated that performing quantum smoothing on a mechanically 
resonant structure when driven by a resonant forcing function gives greater enhancement in precision when 
compared to non-resonant systems. When driven by a Lorentzian process we achieved a simulated improve¬ 
ment in precision of greater than two times better than the equivalent optimal filter, which is consistent 
with theory. We have also experimentally validated the simulation using an experimental testbed. The 
simulations have identified a good parameter regime where greater improvement should be possible in future 
experiments. With future improvements in the system we expect to see further precision enhancements. In 
future work it should be possible to demonstrate further improvement in precision by the incorporation of 
quantum enhancement using a phase squeezed probe beam. These results demonstrate the advantage of 
resonances when performing quantum parameter estimation. This is an initial proof of concept that may 
have applications in areas where mechanical systems are being measured in quantum limited domains. 
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